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We calculate the four-wave mixing (FWM) in a Bose-Einstein condensate system having multiple 
spin wave packets that are initially overlapping in physical space, but have nonvanishing relative 
momentum that cause them to recede from one another. Three receding condensate atom wave 
packets can result in production of a fourth wave packet by the process of FWM due to atom-atom 
interactions. We consider cases where the four final wave packets are composed of 1, 2, 3 and 4 
different internal spin components. FWM with 1- or 2-spin state wave packets is much stronger 
than 3- or 4-spin state FWM, wherein two of the coherent moving BEC wave packets form a 
polarization-grating that rotates the spin projection of the third wave into that of fourth diffracted 
wave (as opposed to the 1- or 2-spin state case where a regular density-grating is responsible for the 
diffraction). Calculations of FWM for 87 Rb and 23 Na condensate systems are presented. 

PACS numbers: PACS Numbers: 03.75.Mn, 03.75.Nt, 03.75.Kk 



I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) 
of dilute alkali gases has created considerable interest 
in the field of nonlinear atom optics. In a BEC, pair- 
wise atomic collisions whose strength is characterized by 
the two-body s-wave scattering length, a, gives rise to a 
nonlinear interaction that can produce four-wave mixing 
(FWM) phenomena. The FWM process allows one to 
study a variety of phenomena ranging from Bose stimu- 
lation P, 0, E|, elastic scattering loss to empty modes 

0, entanglement and correlations ||, and squeezing. 
The first experimental demonstration of FWM of matter- 
waves was reported by Deng, et al. and involved four 
BEC wave packets in identical internal spin states. The 
theory of single spin matter- wave FWM was developed 

in Refs. HHTlll. But FWM is also P ossible for differ " 
ent internal spin states Q , where spin exchange collisions 
may be involved. 

This paper makes predictions about the strength of the 
signal in multiple spin FWM experiments using a mean 
field theory with arbitrary internal spin BEC matter- 
wave packets, including both elastic and inelastic loss 
processes. FWM can occur in BEC systems containing 

1, 2, 3 or 4 different internal spin components. With 
1- or 2-spin states, the process of FWM is analogous to 
Bragg diffraction of matter-waves off a density-grating 
formed by two of the moving BEC wave packets, from 



which the third wave packet can scatter to produce a 
fourth wave. We show that the 3- or 4-spin state case is 
phenomenologically different from the 1- or 2-spin state 
case. In particular, the nonlinear coupling strength of 
the former depends on differences of scattering lengths, 
greatly reducing the population of the fourth wave gen- 
crated by FWM. In this latter case the grating is no 
longer a density-grating but a spin-density grating, or 
polarization-grating, and the diffraction process rotates 
the spin projection of the third wave packet into that of 
a fourth diffracted wave packet. 

This paper is organized as follows. Sec. II provides 
a general description of the process of FWM of multi- 
spin Bosc Einstein Condensates including the Bragg out- 
put coupling technique for generating high momentum 
wave packets. We review general multi-spin wave packet 
coupled equations of motion within the Slowly Varying 
Envelope Approximation (SVEA) and formulate phase 
matching conditions for multi-component spin systems. 
Sec. Ill describes calculations for specific cases of 2-, 3- 
and 4-component cases. Finally, Sec. IV gives a short 
summary and conclusion. 



II. THEORY 



Creating a Moving BEC Wave Packet: 
Conditions 



Initial 
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Fig-fflschematically shows the process of creating mov- 
ing BEC wave packets via Bragg scattering wherein two 
sets of Bragg laser pulses create two moving daughter 
wave packets from an original condensate wave packet. 
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Two laser pulses of central frequency uj and u> + 5 and 
wave vector k w and k w+( 5 incident on a gas of Bosc- 
condcnscd atoms impart a well defined momentum "kick" 
to the atoms. We assume that the frequency u> is close to 
resonance with an atomic transition, and the detuning 
frequency 8 is very small compared with lu. The fre- 
quency u) is chosen to be close to an allowed atomic tran- 
sition (detuned by GHz), and S for the pulse of central 
frequency oj + 5 is chosen to be in the kHz range. The an- 
gle between the propagation directions of the light pulses 
is 9 (0 = 180° corresponds to the counter-propagating 
pulses). We consider a parent condensate with atoms 
in a single internal spin state \Fi,Mi). A set of optical 
light pulses incident on the parent condensate with av- 
erage momentum (Po) = can Bragg scatter atoms via 
Raman scattering from the initial wave packet, thereby 
creating a new daughter wave packet with momentum 
(P) = h(k u — k w+< 5). This process is associated with ab- 
sorption of one photon from the first pulse and stimulated 
emission of one photon into the second light pulse. This 
is a resonant transition and occurs only if conservation 
of energy and momentum are satisfied simultaneously. 
If the resulting velocity acquired by the atom is large 
compared to the speed of sound in the condensate, the 
dispersion relation is quadratic, and conservation of en- 
ergy gives hS = P 2 /2m+Ep ft M f — EFi,Mi , where m is the 
atomic mass, and the subscripts / and i indicate the final 
and initial state respectively. Conservation of momentum 
yields the relation P = 2hk sin(6*/2), where k = 2oj/c\ the 
atomic recoil momentum equals the difference between 
the central momenta of two pulses. It is often convenient 
to specify the velocity of the Bragg scattered atoms in 
units of the recoil velocity vr = %k/m (which equals 2.9 
cm/s for 23 Na). The intensity and duration of the op- 
tical pulses determine the number of atoms that receive 
the recoil "kick" by undergoing the Raman process. The 
polarization and detuning of the lasers determine the fi- 
nal internal spin state \Ff,Mf). Two different sets of 
Bragg pulses are used to create two separate daughter 
wave packets. These, and what remains of the initial 
BEC, are the three nascent wave packets that partici- 
pate in the FWM process. 

The duration of the Bragg pulses is taken to be short 
(fis time-scale) compared to the mean-field evolution of 
the BEC (typically hundreds of /is). As a result, it is 
easy to insure conservation of energy and momentum 
in the Raman process that produces Bragg scattering. 
The whole momentum distribution of the initial BEC 
can therefore participate in the Bragg scattering process 
if the Bragg pulses are sufficiently short. We therefore 
use a "copy" approximation for the initial conditions of 
each slowly-varying envelope created by Bragg scatter- 
ing. The three nascent wave packets can be represented 

as $ r (x, 0) = y^p^ (x, t = 0), for r = 1, 2, 3 at time 

t = 0, < &4(x, 0) = 0, where \^o( x , 0) is the spatial com- 
ponent of the initial wave function, N is a total number 
of atoms and N r (0) is the initial number of atoms in the 




t = t ~ 1 ms 



FIG. 1: Bragg scattering at time t=0 creates two moving 
daughter wave packets from the original parent condensate. 
In the center-of-mass frame in position space, shown in figure, 
all three wave packets are moving with momenta of magni- 
tude h\k\ related to the photon recoil energy. The nonlinear 
interaction of the three initially created wave packets gener- 
ates a fourth wave packet moving in the direction satisfying 
phase-matching criteria. At t — 1 ms the four wave packets 
with spin states \Fj,Mj), j = 1, 2, 3, 4, have moved away from 
one another. 

r-th component. Before the Bragg pulses are applied, 
all the atoms are in the same internal spin state. If the 
Bragg pulse sequence triggering the FWM process is as- 
sociated with a change in the internal spin state of the 
atoms, the daughter wave packet has not only a different 
momentum than the parent, but also a different internal 
spin component. It is essential that the Raman detuning 
from the excited hyperfine state is smaller than the hy- 
pcrfinc splitting in order for Bragg pulses to change the 
internal spin state. 



B. SVEA Equations 

Our zero temperature theoretical model for FWM 
involves condensate dynamics described by the Gross- 
Pitaevskii equation (GPE) We start our analysis of 
different multi-spin component FWM by deriving a set 
of coupled GPEs for all the wave packets participating 
in the process. The Hamiltonian of the system in the 
second quantization can be written as: 

H=J2 [ *«(*>*) I T x + ^( X ^))] *a(x,i)dx+ 

+ \ E / /^(x,<)^(x',i)C/ Q/3 , 75 (x-x') 

a,/3,7,<5 

x * 7 (x',t)* 5 (x,t)dxdx' (1) 
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where subscripts a,/3,j,S denote different spin compo- 
nents. T x = ^-V^ is the kinetic energy operator, 
V(x, t) is an external potential imposed on the atoms 
and U(x — x') is the interaction between particles, which 
for the dilute bose gas, is conventionally taken to be of 
contact form, 

U af3 „ s (x - x') = 4?r ^' 7 ^ (x - x') . (2) 

Here a a p n s is the s-wave scattering length for the mul- 
tispin collision process (see below) and m is the atomic 
mass. The operators vE , Q (x, t) satisfy equal time bosonic 
commutation relations 

[v]/ Q (x, t), ^(x', t)} = [tft ( X) i), vj/t ( x ', *)] = (3) 

[*a(x, t), ^(x', t)] = <5(x - x')<W (4) 

The Heiscnbcrg equations of motion for all *ff a (x, t) fields 
can be obtained by taking the commutator with the 
Hamiltonian QJ. We consider the case when the to- 
tal wave function consists of four wave packets moving 
with different central momenta P r = ?ik r ,r = 1, ... ,4. 
Within a mean field approximation we decompose the 
total wave function into separate wave packets, centered 
about momenta P r 

4 

^(x, t) = 2_, ^r( x j exp(ik r ■ x — iu r t)\a r ) . (5) 

r=l 

where the rth wave packet with internal spin a r and 
mean kinetic energy hu> r — h 2 k 2 /2m: $ r (x, i) is the 
slowly varying envelope (SVE) of packet r, S,k r is the 
central momentum of packet r and \a r ) = \F r ,M r ) la- 
bels the internal atomic spin state of the wave packet 
r. 

Upon substituting the SVE expansion for if? in Eq. (0) 
back into the GPE, collecting terms, multiplying by 
the complex conjugate of the appropriate phase factors, 
and neglecting all terms that are not phase matched, 
thereby making the slowly varying envelope approxima- 
tion (SVEA) [3j, we obtain a set of coupled equations 
for the slowly varying envelopes $ r (r, t). The SVEA was 
necessary to restrict the momentum components only to 
those around the central momentum of each of the wave 
packets in the numerical calculation; otherwise, the mo- 
mentum range would have been too large to treat numeri- 
cally. Solving for the SVE allows efficient numerical simu- 
lations and helps in painting a clear picture of the process 
since it separates out explicitly the fast oscillating phase 
factors representing central momentum Kk r . The SVEs 
$ r (x, t) vary in time and space on much longer scales 
than the phases. The multicomponent SVEA equations 
can be written as 

(| + (-v r ).V + l(-|lv 2 + K(x,t))) $r = 



- -A[/ rl . !rr |(f> r | 2 $ r [self phase] 
--N(l + 5 a p) U rs , rs \<S> s \ 2 <S> r [cross phase] 

s (s^r) 

~Nf aP , Sl J2 CV,»**S*«** [4WM source] 

s,t,q (s=£t=jtq^=r) 

(6) 

where a,0,6,~f are the respective spin components of 
wave packets q,r,s,t, and fap.s-y — 2 if all spin com- 
ponents are the same @(see also Sec. IIII Al below^l. = 
1 if there are two or four distinct spin components, 
and = \/2 if there are three distinct spin components. 
We assume here that the scattering lengths that enter 
into the coupling constants in Eq. (j2J) arc obtained from 
symmetrized scattering matrix elements, as described by 
Stoof et al. H3- 

The interaction of the atoms with an external trap- 
ping field, which could be magnetic or optical, results in 
a harmonic oscillator potential V r (x, t) for each spin com- 
ponent. We allow for a time dependence of the potential 
in order to account for the turning off of the trapping 
field in the experiments we model. 

We only consider the case of zero or extremely weak 
magnetic field, such that the magnetic Zeeman shift is 
very small compared to the recoil energy hk 2 /2m. In 
this case, assuming that the kinetic energy associated 
with the Bragg 'kicks' are much larger than the chemical 
potential of the condensate, the momentum and energy 
conservation conditions have the form 

k r + k 9 -k s -k t = 0, (7) 
k 2 r +k 2 q - k 2 s - k 2 = (in CM frame) . (8) 

Each of the indices r,q,s,t may take values between 1 
and 4. If there were a non-negligible Zeeman splitting 
in the energies of different spin-components, then the en- 
ergy conservation condition would need to be modified 
to take into account the energy released in the inelas- 
tic scattering process. In general, one would expect the 
prospects for FWM to be greatly diminshed if the packets 
have large relative velocities. 

Given our assumption of zero magnetic field, Eqs. I|7)l 
and (JSJ) are automatically satisfied in two cases: (a) r = 
q = s = t (all indices arc equal), or (b) r = s ^ q = t 
(two pairs of equal indices). The corresponding terms 
describe what is called in nonlinear optics self and cross 
phase modulation terms respectively. The self and cross 
phase modulation terms do not involve particle exchange 
between different momentum components. They modify 
both amplitude and phase of the wave packet through 
the mean-field interaction. 

Particle exchange between different momentum wave 
packets occurs only when all four indices in the last term 
on the RHS of Eq. © are different. In this case, satis- 
fying the conservation of momentum and energy is not 
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automatic. 4WM can be viewed as a process in which one 
particle is annihilated in each wave packet belonging to 
an initially populated pair of wave packets and simulta- 
neously one particle is created in each of two wave pack- 
ets of another pair, one of which is initially populated 
and the other is initially unpopulated. Hence, 4WM re- 
moves one atom from each of the "pump" wave packets 
(without loosing generality we may call them 1 and 3), 
and places one atom in the "probe" wave packet (we will 
call it wave packet 2) and one atom in the 4WM output 
(wave packet 4). The bosonic stimulation of scattering 
mimics the stimulated emission of photons in an optical 
medium. The phase matching conditions are particularly 
simple in the center-of-mass frame, defined by the con- 
ditions fci = k 3 , k 2 = ki and k 2 = k\. This picture is a 
consequence of the nature of the nonlinear terms in the 
four SVEA equations. 

For the case of multi-component FWM, there is an 
additional degree of freedom that must be included. In 
this case there is a coupling between different internal 
spin components, introducing new selection rules namely: 

F!=F 2 =F 3 = Fi or F 1 = F 2 ^ F 3 = F 4 , (9) 

and additionally 

Mi + M 3 = Mi + M A . (10) 

The initial wave function is obtained from the numer- 
ical solution of the time-dependent GPE using an imag- 
inary timepropagation in the presence of the magnetic 
potential y|. In the FWM experiments we model, after 
turning off the trapping potential, the BEC wavepacket 
is allowed to freely evlolve for a time ts before the sets 
of Bragg pulses are applied to produce three initial wave 
packets with three different momenta. To model this, we 
propagate the initial wave packet (in real time) for a pe- 
riod of time ig to provide the initial condition in Eq. JSJ). 
This free evolution causes a spatially varying phase to de- 
velop across the condensate as it expands in the absence 
of the trapping potential due to the mean-field interac- 
tion. Given the initial condition, the SVEA equations 
can be used to propagate the envelope function of each 
wave packet, using the same numerical method used to 
propagate the ordinary time-dependent GPE. 

C. Nonlinear Coupling Constants and Loss Terms 

The goal of this paper is to estimate the number of 
atoms in the created FWM wave packet in various re- 
alizations of the multi-spin FWM experiments. To ac- 
complish this goal we must first determine the nonlinear 
coupling constants U qrtS t for 23 Na and 87 Rb in the F = 1 
and F = 2 internal spin states that appear in the GPE. 
We performed the calculations of the strength of the non- 
linear coupling constant U rq ^ s t{k t s), which is determined 
by the following two-body scattering process 

\F S , M s ) + \F u M t ) \F q , M q ) + \F r ,M r ) , (11) 



at relative momentum k st = |k s — k t |. The non- 
linear coupling constant is defined by U qrtS t(k s t)/fi = 
(—T qrjSt /k), where S = l + 2iT is the unitary scatter- 
ing matrix. For the elastic collision case, and for small 
values of relative momentum k st , the coupling constant 
is 

» a st = A st - iB st , (12) 

where a s t is the complex s-wavc scatttering length |llj| . 
The total elastic cross section is given by er c i a stic — (1 + 
8 s t)^'K(A 2 st +i? 2 t ), and the total rate constant for inelastic 
collisions is (i^° tal ) st = E qr ( K 2)st-> qr = f B st . The 
elastic and inelastic collisional losses from the moving 
wave packets can be calculated using er c i as ti c and K 2 , as 
described in Ref. Q. We calculate the various A st and 
B st values for 23 Na and 87 Rb collisions using standard 
coupled channels models of threshold scattering. 

Insight into the nature of the four-wave mixing source 
terms can be obtained by using a scattering representa- 
tion that is useful at low magnetic field, where the total 
angular momentum F t = F1+F3 + I = f + lof the 
colliding pair of atoms is conserved in a collision. The T- 
matrix elements we need are for the case when two atoms 
collide in the internal spin levels (Fi,Mi), (F 3 ,M 3 ) and 
with relative angular momentum 1, to and end up in the 
levels (i^, M2), (F4, M4) with relative angular momentum 
£',m'. We assume collision energy is low enough that only 
s-waves contribute, so that I = £' = 0, m = ml = 0. 

To obtain the relevant quantum numbers, first vec- 
tor couple the |i*iMi) and I.F3M3) levels to obtain 
the resultant angular momentum \F\F 3 f,M\ + M3) 
states, then vector couple the relative angular momen- 
tum |£, to) to get the total angular momentum states 
\F x F 3 f,i,F t ,M t = Mi + M 3 + m). The desired sym- 
metrized T-matrix elements 01 are given in terms of 
the indices {FiF 3 f£F t } and are independent of projec- 
tion quantum numbers Mi, M3, to, M t . 

The only significant collisions we need to consider are 
spin-exchange collisions. These are possible only if the 
following selection rules are obeyed: I = to = to', 
/ = /', and Mi + M 3 = M 2 + M 4 ; for s-waves i = 
and f = f = F t . The four- wave mixing source term 
is proportional to a(st — > qr) = —T qr ^ t /k s t. Using the 
transformation of T-matrix elements to the total angular 
momentum basis 0, assuming I = £' = we 

obtain 

aCfi-Mi, .F3M3 F 2 M 2 ,F 4 M 4 ) = 

f 1 + S FlFa \ * / 1 + 5 F2Fi \ * 

\l + Sf 1 F 3 S]\i 1 m s ) \ 1 + $F 2 F4. Sm 2 M 4 / 

x ^{F 2 F 4 f\M 2 Mi, M 2 + M 4 ) 

/ 

x(F 1 F 3 f\M 1 M 3 ,M 1 +M 3 )aV\F 1 F 3 -» F 2 F 4 ), 

(13) 

where ( | ) are Clebsch-Gordan coefficients. Clearly four- 
wave mixing is possible only if the collisions are energy 
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degenerate: F 2 = F\ and F4 = F 3 . 

For 1- and 2-component four wave mixing, where Mi = 
M 3 = M 2 = Mi or M 2 = A/i and M 4 = M 3 respectively, 
the product of Clebsh-Gordon coefficients is positive def- 
inite, and the sum will be of the same magnitude as the 
individual elements. On the other hand, for three- or 
four-component four wave mixing, there will be terms 
with the product of the Clebsh-Gordon coefficients being 
positive and terms with the product being negative. In 
this case, the scattering length that controls the FWM 
source term in Eq. © will involve differences of scatter- 
ing lengths of comparable magnitudes, and will tend to 
be much smaller than for the one-or two-component case. 
For example, 

a(l - 1,11 ^ 10,10) = 

^(a (2) (ll-ll)-« (0) (ll-H)) , (14) 
a(l - 1,2- 1 -> 10,2- 2) = 

^j- (a (3) (12 -» 12) - a (2) (12 -> 12)) . (15) 

Inelastic energy-releasing exchange collisions arc only 
possible for the cases{i^F 3 /} = {122}, {220} and {222} 
for 23 Na and 87 Rb. For 87 Rb, these loss collisions are 
anomously small 0,0,0 and we can set (-F1-F3) = 
0. 

Table H] gives our calculated scattering lengths, based 
on coupled channels models of threshold collisions, that 
should be used in the U qriS t coupling constant that gives 
the source term for FWM in Eq. i.e., U qr<s t = 

Antaqr sth /m. We include several cases for 23 Na and 
87 Rb involving one, two, three, and four spins. The three 
and four spin cases have significantly larger coupling con- 
stants for 23 Na than for 87 Rb, because scattering length 
differences in (|15|) tend to be much larger for the for- 
mer case. The imaginary part of the complex scattering 
length for the source term is small compared to the real 
part, even for 23 Na collisions, and is not shown. 



III. RESULT OF CALCULATIONS 

A. One- And Two-Spin Component FWM 

In the numerical simulations we model the condensate 
comprised of magnetically trapped atoms without a dis- 
cernible non-condensed fraction. We assume trap fre- 
quencies in the x, y and z directions of 84, 59 and 42 
Hz, as in Ref.Q. After the magnetic trap is switched off, 
the condensate expands freely during a delay time tg. 
Then the sequence of two sets of Bragg pulses creates 
two moving wave packets 2 and 3 in the lab frame, in 
addition to the initial stationary wave packet 1. One can 
change the momentum of the wave packets by changing 
the angle between the laser beams used for the Bragg 
out-coupling or by changing the laser frequencies. Our 



simulations neglect the detailed dynamics during the ap- 
plication of the Bragg pulses; instead, we assume that 
after the Bragg pulses are applied, each wave packet is a 
copy of the parent condensate wave function. We shall 
consider the case where the initial three wave packets 
contain an equal number of atoms (except for the case 
corresponding to the original NIST experiments where 
the ratio is 3:7:3). 

In the center-of-mass reference frame, all the wave 
packets (including the new FWM wave packet) move with 
the same velocity. In the one and two-spin component 
cases the thresholds of the incoming and outgoing colli- 
sion pairs coincide, i.e., F q + F r = F s + F t . Moreover, 
Mp is a good quantum number in the two-body collision 
requiring M q + M r = M s + M t for all cases. The three 
initial wave packets can each be in an arbitrary internal 
spin state of the alkali atom (for 23 Na or 87 Rb this is 
either F = 1 of F = 2). The FWM process preferentially 
populates the fourth wave packet in a spin state that 
satisfies the energy and angular momentum projection 
constraints. 

Of all the possible combinations of spin states we can 
start with, we find two distinct classes of combinations: 
(1) at least two of the three initial wave packets are in 
the same spin state, (2) all three wave packets are in 
different spin states. In both cases, the overlap of the 
coherent moving BEC packets, 1 and 2, form a grat- 
ing (either a density-grating or a spin-density grating) 
that 3 diffracts off, producing a new wave packet 4. In 
the single-spin case the number of atoms in the fourth 
wave packet is four times larger than in the 2-spin case, 
(iV^.spin ~ 4(A r 4)2- S pin even if the scattering lengths 
and initial number of atoms for these cases were exactly 
identical. The grating picture that explains the factor of 
four is as follows. For the 1-spin state case, wave packets 
1 and 2 form a grating and wave packet 3 can scatter 
off the grating to produce atoms that are in wave packet 
4 (see Fig. ^| , and wave packets 2 and 3 form a grating 
and wave packet 1 can scatter off the grating to produce 
atoms that are in wavepacket 4. These two amplitudes 
add up coherently so the probability is four times the 
probability that would be obtained were there only one 
amplitude for producing the 4th wave. For the two-spin 
state case, if packet 1 (or 3) is the different spin state, 
it will scatter off the grating formed by 2 and 3 (1 and 
2) and form the same spin state in a wave packet with 
momentum h^. The gratings are dynamical ones that 
change in time as atoms are removed from wave packets 
1 and 3 and placed in wave packets 2 and 4. 

The s-wave scattering lengths were used to form the 
nonlinear coupling constants, U qr ^ s t = Ana qr ^ st Ti 2 /m, and 
in turn, these were used in our numerical simulations of 
FWM. Fig. |21 shows the fraction of atoms in the newly 
generated wave packet, = N4/N, versus the total num- 
ber of atoms in the parent condensate N for a number 
of different FWM processes. We assumed a free expan- 
sion time (i.e., the delay time) of f# = 600 fis, as in 
the NIST experiment |lj. The strongest FWM convcr- 
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sion is obtained for |1, + 12, 0} 3 — > |1, -1) 2 + |2,0) 4 
collisions of 87 Rb atoms, since this case has the largest s- 
wave scattering length in the source term. The saturation 
(and even decrease in the 87 Rb case) of the FWM fraction 
= TV4/TV with increasing TV is clear from Fig. The 
origin of this saturation is discussed in Ref. 0] and is due 
to the physical separation of the interacting wave pack- 
ets, back conversion from the newly formed wave packet 
4, and elastic and inelastic scattering loss processes. 

For reference, we included in Fig. [21 the original data 
from the NIST experiment []]. Comparison of the ex- 
perimental results for 23 Na |1, — 1) with the theoretical 
results show that at the highest values of TV, the calcu- 
lated number of atoms in the FWM wave packet is higher 
than in experiment. It is possible that stimulated elastic 
scattering loss may have to be taken into account at large 
TV. 
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FIG. 2: Fractional FWM output versus total number of atoms 
N, /4 = TV4/TV, in cases when one- or two-spin components 
are present. The simulations were carried out allowing the 
condensate to expand for 600 /is before applying the Bragg 
pulses, as in the NIST experiment QJ. The solid circles are the 
experimental results obtained for 23 Na |1, — 1). In the two- 
spin component simulations the ratio of atoms in the three 
initially populated wave packets are 1:1:1, whereas the one- 
spin component simulation used the experimental ratios of 
7:3:7. 



The scaling of the fraction of atoms output into the 
FWM wave packet with N 2 (t )/N for TV = TVi(t ) + 
TV 2 (t ) + TV3O0) = TV = 1.5 x 10 6 is shown in Fig. for 
the case of 1 1, -l)i + 1 1, 0) 3 — >|1,-1) 2 + |1,0) 4 collisions 
of 23 Na atoms. Note that on the left part of the figure 
TVi(to) = TV 3 (i ) > N 2 {t ), while on the right part of 
the figure TVi(i ) = TV 2 (t ) > TV 3 (t ). In Ref. ||, we 
discussed how the fraction of atoms in the FWM wave 
packet would scale with the initial number of atoms in 
the three nascent wave packets. There we showed how 
a simple argument predicts that the total FWM output 



fraction is given by 

TV4 (t coi ) 



N 



( t C ol 
\tNL 



(16) 



where /, = TVj(fo)/TV, t co i is the time for the wave pack- 
ets to separate, and t^L is the characteristic nonlinear 
time scale, t^^ = h//i where /i is the chemical poten- 
tial. This is an upper bound on the 4WM output, since 
the mutual interaction of the packets due to the self- 
and cross-phase modulation terms (the self- and cross- 
interaction energy terms), and the elastic and inelastic 
loss processes are not included in the estimate. The curve 
labeled TViTV 2 TV 3 /TV 3 in Figgis obtained using Eq. i|16[l . 
It follows the calculated results rather well. Clearly, for 
cases where TV < 1.5 x 10 6 , one expects that the simple 
argument will Drovide as scood an estimate. 

23 Na (11,-1^ + |1,0> 3 ->|1,-1> 2 + |1,0> 4 ) 



n = N^y+Nyy+isyy+N^y = 1.5x10° 

— N^Nj/N 3 scaling 




FIG. 3: Fractional FWM output versus N 2 (t )/N. The solid 
curve is the result of using the simple model in Eq. I|16|l and 
the triangles are the result of calculations. The trap is the 
same as used in Fig. |2] The Bragg pulses are applied 300 /is 
after the trapping potential is turned off. 

The dependence of the fraction of atoms in the FWM 
wave packet on the velocity of the condensate wave pack- 
ets is shown in Fig. 0] for the case of 23 Na |1, — l)i + 
1 1, 0)3 — > |1,— 1)2 + |1,0)4 condensate collisions where 
the free-expansion time tg between trap off and the 
Bragg pulses is 300 fxs. The velocities of the wave pack- 
ets in the center-of-mass frame are indicated in Fig. 0] in 
units of the recoil velocity vr (recall that the recoil ve- 
locity is defined as vr = Kk/m which equals 2.9 cm/s for 
23 Na). The relative velocities of wave packets 1 and 3 or 
2 and 4 are twice these velocities. As the relative veloc- 
ity increases, the number of atoms in the FWM packet 
decreases, since the duration of the wave packet overlap 
essential to FWM decreases. For example, for the curve 
labeled I.Ovr in Fig. 0] it takes about 500 /us for the cen- 
ters of the moving wave packets to separate by the mean 
Thomas-Fermi diameter of the initial condensate. The 
saturation of TV4/TV is evident. For very slow relative 
velocities, TV4/TV actually begins to decrease with time 
due to back-conversion via FWM, as first described in 
Ref. 0. 
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Na ( 11,-1^ + I1,0> 3 - ll,-l> 2 + ll,0> 4 ) 
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FIG. 4: N 4 (t)/N versus t for N = 1.5 x 10 6 atoms for various 
velocities hk/m of the wave packets. The trap is the same as 
used in Fig. [5] The Bragg pulses are applied 300 /is after the 
trapping potential is turned off. 



A substantial loss of the condensate atoms due to elastic 
and inelastic scattering collisions occurs, but nevertheless 
the newly created FWM wave packet contains 120,000 
atoms after full separation. This number of atoms can 
be easily detected. In fact, for all the cases shown in 
Fig. the generated FWM signal is strong enough to be 
detected in real experiments when N > 1.0 x 10 6 atoms. 




1500 



The dependence of the fraction of atoms in the FWM 
wave packet on the free-expansion time tE between turn- 
ing off the harmonic potential and applying the Bragg 
pulses is shown in Fig. [S] for the case of 23 Na |1, — l)i + 
1 1,0)3 — ► 1 1, — 1)2 + |1,0)4 condensate collisions when 
the wave packets velocity is vr = 2.9 cm/s. During 
the free-expansion of the condesate, a spatially varying 
phase develops across the parent condensate, and this 
phase dclctcriously affects the phase matching required 
for FWM. This has already been discussed in the single- 
component studies that were published previously 

23 Na ( ll,-l>j + ll,0> 3 -»- Na ll,-l> 2 + H,0> 4 ) 



0.06 



0.04 - 



0.02 - 



N = 1.5x10° 



T = |xs 




T = 300 fis 



T F = 600 us 



600 



800 



FIG. 5: N 4 (t)/N versus t for N = 1.5 x 10 6 atoms for various 
free-expansion times tE- The trap frequencies are the same 
as used in Fig. [5] 

Figure shows the total number of atoms in all the 
wave packets and the number of atoms in each of the 
wave packets during the FWM half-collision versus time 
for the 87 Rb |l,-l)i + |2,0) 3 — ► [1, — 1> 2 + |2,0) 4 case. 



FIG. 6: N(t), Ni(t), N 2 (t), N 3 (t), and N 4 (t) versus t for 
87 Rb |l,-l)i + |2,0> 3 — > |1,-1)2 + |2, 0)4 collisions with 
TV = 1.0 x 10 6 atoms. The trap frequencies are the same as 
used in Fig. |2] The Bragg pulses are applied 600 /is after the 
trapping potential is turned off. 



B. Three- and Four-Spin Component FWM 

FWM is possible with any combination of arbitrary in- 
ternal spin states, provided Me is conserved in the elastic 
two-body scattering process. FWM processes cannot be 
interpreted as Bragg diffraction off density-gratings in 3- 
or 4-spin FWM. Instead, a "polarization-grating" scat- 
ters the atoms, as described above. For example, the 
four-mixing output (packet 4) in the process 



11,-1) 



|1,+1) 



|1,0) 2 + |1,0) 4 



(17) 



can be thought of as a rotation of the spin "polarization" 
of wave packet 3 due to scattering off the spin density 
grating formed by the overlap wave packets 1 and 2. 

Fig. □ shows Ni/N versus N for 23 Na |l,-l)i + 
|1,+1) 3 -> 1 1, 0)2 + |1,0} 4 and for 23 Na |l,-l)i + 
|2, —1)3 — ► |1,0)2 + |1, ^2)4 condensate collisions. Fig. El 
shows the time-dependent population fractions for the 
latter 4-component case for TV = 2.0 x 10 6 atoms. The s- 
wave scattering length in the source term for the former 
process is a = 0.14 nm, whereas for the latter process, 
a = 0.59 nm. Taking an angle of 20 degrees between 
the Bragg laser pulses ensures a relatively low velocity of 
hk/m = 1 cm/sec for the moving packets. This allows 
the packets to remain overlapped long enough to generate 
an observable fourth wave. 
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FIG. 7: Fractional FWM output versus total number atoms 
JV calculated upon allowing the condensate to expand for 
600 fis before applying the Bragg pulses for 23 Na |1, — l)i + 
|1,+1) 3 |1,0> 2 + |1,0> 4 and for 23 Na |1, -l) x + |2, -1) 3 
|1, 0)2 + |2, —2)4. The initial ratio of atoms in the three pop- 
ulated wave packets are 1:1:1. The trap frequencies are the 
sa □ 



1 1 1 

M.o> 2 




^^_il'- 1> 1 








|2,-2> 4 

i , i 



1000 2000 3000 



FIG. 8: Ni(t), N 2 (t), JV 3 (t), and N 4 (t) versus t for JV = 2x 10 6 
atoms and 23 Na |1, -l)j + |2, -1) 3 -> |1, 0> 2 + |2, -2) 4 . The 
trap frequencies are the same as used in Fig. [5] The Bragg 
pulses are applied 600 [is after the trapping potential is turned 
off. 



IV. SUMMARY AND CONCLUSIONS 

We have developed a general theory for describing four- 
wave mixing (FWM) of matter-waves in arbitrary inter- 



nal spin states within the context of a mean-field theory 
using the Gross-Pitaevskii equation (GPE). The slowly 
varying envelope approximation is used to write sepa- 
rate equations for each of the condensate wave packets. 
The atom-atom interactions of the BEC atoms result in 
both the mean-field interaction terms, analogous to the 
self and cross phase modulation terms in nonlinear op- 
tics, and the FWM source terms in the GPE. The the- 
ory also incorporates elastic and inelastic scattering loss 
processes. These processes take atoms out of the conden- 
sate wave packets and therefore reduce the FWM. FWM 
with 1- or 2-spin states is analogous to Bragg diffrac- 
tion of matter-waves off a density-grating formed by the 
moving BEC wave packets. The 2-spin state FWM is 
generally smaller than the 1-spin state case by a factor 
of about four due to the coherent addition of two am- 
plitudes for scattering off the density gratings formed in 
the 1-spin case, whereas only one scattering amplitude 
occurs for the 2-spin state case. FWM with 3- or 4-spin 
state are generally much weaker processes; in the 3- or 4- 
spin state cases, the coherent moving BEC wave packets 
form a polarization-grating (a spin-density grating) that 
rotates the spin projection of the diffracted wave. 

Calculations of multicomponent FWM for 87 Rb and 
23 Na condensate systems have been presented. In these 
calculations the 3- and 4-spin state FWM output signals 
are lower by roughly an order of magnitude than for 1- 
and 2-spin state FWM cases. The reduction is due to the 
much smaller source term for FWM in the 3- and 4-spin 
state cases, since the coupling strength involves differ- 
ences of scattering lengths (see Eq. 1)15(1) of comparable 
magnitudes. 
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